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ABSTRACT 

Magnetic fields are dragged in from tlie interstellar medium during the gravi- 
tational collapse that forms star/disk systems. Consideration of mean field mag- 
netohydrodynamics (MHD) in these disks shows that magnetic effects produce 
subkeplerian rotation curves and truncate the inner disk. This letter explores 
the ramifications of these predicted disk properties for the migration of extraso- 
lar planets. Subkeplerian flow in gaseous disks drives a new migration mechanism 
for embedded planets and modifies the gap opening processes for larger planets. 
This subkeplerian migration mechanism dominates over Type I migration for suf- 
ficiently small planets (mp < IM®) and/or close orbits (r < lAU). Although 
the inclusion of subkeplerian torques shortens the total migration time by only 
a moderate amount, the mass accreted by migrating planetary cores is signifi- 
cantly reduced. Truncation of the inner disk edge (for typical system parameters) 
naturally explains final planetary orbits with periods P ~ 4 days. Planets with 
shorter periods P ~ 2 days can be explained by migration during FU-Ori out- 
bursts, when the mass accretion rate is high and the disk edge moves inward. 
Finally, the midplane density is greatly increased at the inner truncation point of 
the disk (the X-point); this enhancement, in conjunction with continuing flow of 
gas and solids through the region, supports the in situ formation of giant planets. 



Subject headings: planetary systems: formation — solar system: formation — 
stars: formation — turbulence — MHD 
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1. INTRODUCTION 

The population o f observed extraso lar planets exceeds 350, with many objects found 



in short-period orbits (lUdry et al. 20071 ). Tremendous progress has been made in our the- 
oretical understanding of the migration mechanisms that allow planets to move (usually 
inward) as they form (Papaloizou & Terquem 2006, hereafter PT06). Parallel advan ces have 



been made in our understanding of star formation processes (IReipurth et al. 20071 ). which 
produce the circumstellar disks that provide birth places for planets. However, these latter 
developments have not been fully integrated into theories of planet formation/migration. 
This letter explores the implications of the predicted disk properties for planet migration. 

During the collapse phase of star formation, magnetic fields are dragged inward and 
retained by the nascent disk. These magnetic fields, along with those generated by the star, 
can drive powerful outflows (Blandford fc Payne 1982, Shu et al. 1994), provide channels for 



accretion of disk material onto the star (IBouvier et al. 20071 ). and produce a back reaction 



on the structure (Shu et al. 2007, hereafter S07). In particular, circumstellar disks are 
predicted to display subkeplerian rotation curves. We parameterize the rotation speed of 
the gas through the ansatz Qg^s = f^, where = GM^/r^. Recent MHD calculations 
indicate that / ~ 0.66 for disks around T Tauri stars and / ~ 0.39 for disks undergoing 
FU-Orionis outbursts (S07). The relative velocity between the gas and orbiting planets is 
thus a substantial fraction of the Keplerian rotation speed, and / 7^ 1 over much of the disk. 

Planetary migration is altered when the circumstellar gas moves at subkeplerian speeds. 
In Keplerian disks, embedded planets drive wakes in to the disk, and asymmetric g ravitational 



interactions of these wakes drive planets inward (iGoldreich fc Tremaine 19791 ). a process 
called Type I migration (Ward 1997ab). In subkeplerian disks, embedded planets orbit with 
Keplerian speeds and experience a headwind from the gaseous disk (which moves slower). 
This velocity mismatch results in energy loss from the orbit and inward migration. This paper 
explores this mechanism (Section 2), including an assessment of the conditions necessary for it 
to compete with Type I migration. A weaker version of this effect occurs in nearly Keplerian 
disks, where 1 — / = 0{H'^/r'^) ~ 1/400, where H is the disk scale height. The headwind 
produced by this small departure from Keple rian flow affects dust gr ains as they coagulate 



into larger bodies, but does not affect planets (IWeidenschilhng 1977bl ). The departures from 



Keplerian speeds due to magnetic effects are larger by two orders of magnitude. 

The second feature of magnetic disks is that their inner edge is truncated. Disk material 
accretes onto the star by moving along magnetic field lines in a "funnel flow" . In one theory, 
the disk truncation point, the start of the funnel flow, and the origin point for outflows 
occur at a location known as the X-point (Cai et al. 2008, hereafter C08). This "point" 
has a finite extent, comparable to the disk scale height. In competing theories, the wind 
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originates from an extended disk region (jPudritz &: Norman 19831 ). but the disk retains an 
inner truncation point; disk wind models also predict subkeplerian flow, which might inhibit 
wind launching (Shu et al. 2008). Because the disk does not extend to the star, migration 
is naturally halted at the truncation point (r ~ 0.05AU), as recognized (ILin et al. 19961 ) 
soon after the discovery of extrasolar planets (Mayor & Queloz 1995, Marcy & Butler 1996). 
This paper argues that observed ultra-short-period planets (P ~ 1 — 2day, r ~ 0.02AU) 
can be understood through disk truncation at the X-point combined with FU-Ori outbursts 
(Section 3). 



2. SUBKEPLERIAN MIGRATION 



This section presents a migration mechanism driven by subkeplerian flow in magnetically 
controlled disks. To illustrate the properties of this mechanism, we consider power-law disks, 
where the surface density and temperature distributions take the form E(r) = Si(ri/r)P and 
T(r) = Ti{ri/rY. We take ri = lAU, so the coefficients Si and Ti correspond to val- 
ues at lAU. The power-l aw index for the surface density lies in the range 1/2 < p < 2 
( jCassen fc Moosman 198ll ). For the Minimum Mas s Solar Nebula (MMSN), the surface 
density constant Si ^ 4500 g/cm^ (IKuchner 20041 ). The power-law index for the tem- 



perature profile lies in the r ange 1/2 < q < 3/ 4. Viscous accretion disks (jPringle 198ll ) 
and fiat reprocessing disks (lAdams fc Shu 19861 ) hav e q ~ 3/4; flared repr ocessing disks 
( IChiang fc Goldreich 19971 ) and the early solar nebula (IWeidenschilling 1977al ) have q ~ 1/2. 
Planetary orbits are taken to be circular, with orbital angular momentum J = mp{GM^rY/'^ . 
The disk scale height H = Cg/fi, where Cg is the sound speed, which is determined by the 
disk temperature profile. The radial dependence of H/r is thus H/r = {H /r)i(r /ri)^^"''^^'^ , 
where {H/r)i ^ 1/20 (S07). 



2.1. Torques 

For circular orbits, the torque Tx exerted on the planet by gas in the subkeplerian disk 
takes the form 

Tx = Cz)7^i?|rpgasWrel^ (1) 

where r is the orbital radius, v^ei is the relative velocity between the gas and planet, pgas is the 
gas dens ity, and Rp is the plan etary radius. The dimensionless drag coefficient Cd ~ 1 for 



planets (IWeidenschilling 1977bl ). but can be enhanced by gravitational focusing (Ohtsuki et 
al. 1988). The gas density Pgas ~ T,/2H. Because the gaseous disk orbits with subkeplerian 
speeds f2gas = f^, the relative velocity takes the form fj-ei = (1 — f){GM^/rY^'^, and the 
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torque becomes 



Tx = lCn{l - ffGM^Rp' (I ) . (2) 



The time evolution of the semi-major axis obeys the differential equation 

Idr _ 2Tx 
r dt mpVlr"^ ' 

where we assume that the orbit remains circular. The migration timescale is given by 



(3) 



For typical parameters, M* = I.OMq, mp = l.OM®, Rp = l.Oi?®, Si = 4500 g/cm 
f = 0.66, and {H/r)i = 0.05, the timescale to ^ 70,000yr. 



For comparison. Type I torque strengths are given by 

2 



T,^C,(^^).^rHrny-{L)\ (5) 

where the dimensionless constant Cj ~ (1.364 + 0.541p)/7r ^ 0.606 (Ward 1997a, Tanaka 
2002). The ratio of the torque Tx to the Type I torque is given by 

2 / n \ 2 



Tx 
Ti 



Cnjl-ff 
2Cr 



MA (RA (H) ocM^m-//\-^'^^y\ (6) 



rrip ) \ T ) \r 

For an Earth-like planet (mp = IM^, Rp = 1-R©) orbiting a solar mass star, Cd = Cj, 
and H/r = 1/20, the two torques have almost the same magnitude at r = lAU when 
(1 — /) = 0.34. Since the torque ratio grows with decreasing semi-major axis, subkeplerian 
torques Tx are stronger for r < lAU. The ratio decreases with the planetary mass, so that 
larger planets are more affected by Type I migration. This torque ratio also depends on the 
stellar mass, so that subkeplerian torques are more important for systems associated with 
more massive stars. Figure [1] shows the regions of parameter space {r,mp, M^) for which 
subkeplerian migration dominates over Type I migration. Subkeplerian torques are relatively 
stronger in the inner disk, where magnetic effects (/ 7^ 1) are most likely to arise. 

Note that this treatment is approximate: The T/ expression neglects magnetic effects, 
which can be important in these systems (Terquem 2003). In general, T/ should be a function 
of /. Moreover, this Type I torque is derived for isothermal disks, whereas X-wind disks have 
temperature gradients, which can affect migration (Paardekooper & Mellema 2006). Finally, 
the two torques (Tj,Tx) are not fully independent — both mechanisms produce disk wakes 
that can interfere with each other. The comparison shown in Figure [1] is thus subject to 
future corrections. 
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Fig. 1. — Regions of {r,mp) parameter space for which migration from subkeplerian flow 
(with / = 0.66) dominates Type I migration. Planetary mass (M®) is plotted versus orbital 
radius (AU). Subkeplerian migration dominates in the region below the lines, which are 
drawn for stellar masses M^/Mq — 2, 1, 0.5, 0.25 (top to bottom). 
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2.2. Gaps 

In standard treatments of gap formation, parcels of gas flowing past the planet are 
scattered by its gravitational potential, and deflected through the angle 6 = {2Gmp/bvrei^), 
where b is the impact parameter of the encounter (PT06). This deflection changes the parcel's 
specific angular momentum by the increment 

where the constant K includes corrections arising from the rotating reference frame (Gol- 
dreich & Tremaine 1980). The disk angular momentum changes at a rate determined by the 
integral of AJ. 

In the Keplerian case, the planet moves faster (slower) than disk material on its out- 
side (inside). The planet thus transfers angular momentum to the outer disk and negative 
angular momentum to the inner disk, and opens a gap from both directions. However, in 
a subkeplerian disk, the planet moves faster than the disk material on both sides of the 
potential gap; this complication modifies the gap formation process. 

Although supersonic flow {v^^i > Cs) alters this picture, we start by estimating the 
angular momentum transferred by "scattering" gas parcels and compare the result to that 
transferred by viscosity. Consider the disk exterior to the planetary orbit. The planet moves 
faster than the gas, and angular momentum is transferred outward at the rate 

where the integral's support is concentrated (by assumption) near the minimum impact 
parameter &min- This expression neglects disk shear terms, which must be included to recover 
the standard result in the limit / ^ 1 (compare equation [8] with PT06). 

Gap formation requires that the rate of angular momentum transfer to the disk from 
"scattering" exceeds the rate at which angular momentum is transferred by viscosity. This 
condition requires 

^ > Svrz/Sr^fi, (9) 

where the viscosity u is parameterized as z/ = {2/3)aCsH , and the viscosity coefficient a is 
dimensionless. 

Gap formation requires a second condition: The planet's sphere of influence must be 
large enough, e.g., the Roche radius Rr = r(mp/3M^,)^/^ must exceed the scale height H. 
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The requirement Rr > H places another constraint on the planet mass required for gap 
opening: 



This condition allows the flow to experience n onlinear deviations (I Ward 1997a[) . including 



shock formation in the vicinity of the planet (iKorycansky fc Papaloizou 19961 ). Since this 
requirement could be modified for high Mach numbers, equation (ITOll should be considered 
as approximate. This constraint allows us to use 6min = Rr oc r to evaluate the angular 
momentum transfer rate in equation ([8]), and solve equation ([9]) to find the planet mass 
required for gap opening: 

[wj - — — "UJ • ^^^^ 

Gaps form when both constraints (ITOl) and (fTTj) are satisfied (see also Crida et al. 2006). 
For a f=^ 10~^, / = 0.66, and H/r = 1/20, the nonlinear constraint (fTOj) and the viscous 
constraint (fTT!) are comparable, and the mass required for gap clearing is ~ lOOM^ (about 
one Saturn mass). Because the planet moves faster than the gas on both sides of the gap, 
the planet readily transfers angular momentum to the outer gap edge, but tends to draw 
disk material outward from the inner edge. When the planet has sufficient mass (equation 
[To]), this material shocks against the wake of the planet. The planet will thus be located 
closer to the inner gap edge. 



2.3. Orbital Evolution 

To consider orbital evolution of the planet, we define a dimensionless angular momentum 
variable ^ = J/Ji, where J = mpVLr"^ and Ji is the orbital angular momentum at r = lAU. 
We also define reduced torque parameters through the ansatz F/ = Ti{Ji)/Ji and Fx = 
Tx{Ji)/ Ji- The equation of motion then becomes 

^ = - [TrC^ + TxC'] , (12) 

where the indices a = 2{p — q) and b = 2p + 3 — q. For typical values p = 3/4 and q = 1/2, 
the torque indices become a = 1/2 and 6 = 4. Since the two torques are nearly equal for 
Earth-mass planets at lAU, Tx ^Tj ^ 10(Myr)-i = l/(0.1Myr). Keep in mind that both 
torque expressions are approximate (see Section 2.1). 

For subkeplerian torques acting alone, the time required to migrate from a starting 
location = Jo/Ji = (ro/lAU)"*^/^ to the central star is given by Txtx = ^o^^/ib + 1). 
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Similarly, the time required for Type I torques to move the planet from ,^0 to the star is 
given by F/t/ = C,Q~^^/{a + 1). Since we expect b ^ a, subkeplerian migration is much 
faster than Type I migration for planets starting inside ~ lAU. For planets starting at 
larger radii, however. Type I torques dominate until the planet reaches r ~ lAU, where 
subkeplerian torques take over. 

The above results assume that the planetary core mass remains constant. However, 
the core mass must increase with time for giant planets to form, and the torque strengths 
change as the mass increases. To illustrate this behavior, we use a simplified treatment of 
planetesimal accretion where the planetary core grows according to 

^ = Q^TrRpRnEzn, (13) 

where oc S is the surface density in solids (Papaloizou & Terquem 1999, Laughlin et al. 
2004a). Giant planet formation typically requires ?a 130 gcm~^ at r = lAU, somewhat 
larger than the MMSN value (Lissauer & Stevenson 2007). For a given density of planetary 
cores, pc ~ 3.2 gcm~^, the planetary radius Rp is determined by the mass. After defining 
a scaled radius t] = Rp/Rq), we write the core accretion process in dimensionless form, 

^ = ir^'^^'\ (14) 

where 7 ~ 20(Myr)~^. We can combine planetary growth with orbital evolution. As the 
planet grows according to equation (IT^ . the torques vary with planetary radius, and the 
equation of motion for ^ = (r/lAU)^/^ becomes 



dt 



where f/ = Tj{ri = 1) and fx = ^x{v = !)• 

When migration reaches r < lAU, subkeplerian torques become larger than Type I 
torques and planetary migration takes place more rapidly. If migration starts at larger radii, 
however, most of the evolution time is spent at larger radii where Type I migration torques 
prevail. Consequently, subkeplerian torques don't substantially change the total migration 
time for starting radii beyond r ~ 2AU. However, the mass accreted during the migration 
epoch is significantly altered due to acceleration of the latter stages. These trends are 
illustrated by Figure [21 which shows migration times and final masses for planets starting 
at varying locations tq with fixed initial mass mp = IM^. Here, the equations of motion 
fll4fl5l) are integrated using the standard parameter values quoted above. 

For completeness, we note that subkeplerian flow o ccurs only in in agnetically active disk 



regions. The flow becomes Keplerian in "dead zones" (iGammie 19961 ). where this migration 
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0.1 1 10 



ro (AU) 

Fig. 2. — Migration time and final core mass versus starting radius. Solid curve shows 
migration time tM (Myr, left axis) including subkeplerian and Type I torques. Dotted curve 
shows Im for Type I torques only. Dashed curves show final core mass mp (M®, right axis) 
for migration with both torques (bottom) and Type I torques only (top). 



- 10 - 



mechanism will not operate. However, ordinary Type I migration continues to act within 
dead zones, so that migration occurs at the standard Type I rate. 



3. DISK TRUNCATION 

This section considers the implications of disk truncation. In X-wind theory, this inner 
cutoff radius is given by 

1/7 

rx= , (16) 



where /i* is the magnetic dipole moment of the star, M^, is the stellar mass, Md is the mass 
accretion rate, and $dx ~ 1 is a dimensionless parameter (Ostriker & Shu 1995, S07). Be- 
cause fx marks the inner disk boundary, planets cannot migrate inside the X-point through 
standard migration mechanisms involving disk torques because there is no disk material. 
The X-point location rx ~ 0.05AU for typical parameters (COS), and is coincident with the 
semi-major axes of "hot Jupiters" . Disk truncation through magnetic effects thus provides 
a natural mechanism for ending migration and explains the observed planets in ~ 4-day 
orbits. However, extrasolar planets have recently been discovered with significantly smaller 
semi-major axe^j, r ^ 0.02AU. Although some variation in rx is expected, the presence of 
extrasolar planets in these extremely tight orbits represents some tension with expectations 
from X-wind theory, and suggests the possibility of an additional mechanism. 

The presence of these short-period planets can be understood as follows: The location 
of the X-point depends on the mass accretion rate t hrough the disk. During FU-Ori outburst 
events, the mass accretion rate is greatly enhanced (jHartmann 2001 ). and the X-point moves 



inward (SOT). Equation (fT6l) indicates that must increase by a factor ~ 25 to move the 
X-point to rx ~ 0.02AU as required. Since outburst phases have Md f» lO~^M0yr~^, planets 
can (in principle) migrate inward during these FU-Ori phases. 

This simple picture has two complications. First, FU-Ori outbursts are relatively short- 
lived. The time-averaged mass accretion rate is Md ~ 10^"^ MqJt^^ . During outbursts, the 
Md is ~ 100 times larger, so that the outburst duty cycle is correspondingly shorter, and 
each event lasts only 100 — 1000 years. During the outburst phase, the disk extends inward 
to the desired radii, and migration can take place. However, the planets will not necessarily 
have enough time to migrate inward far enough during the (short) outburst phases. This 
issue is made more urgent by the timing of the events: The planets take a relatively long time 



^Schneider, J. 2009, Extrasolar Planets Encyclopedia, 



http: / /exoplanet.eu/ catalog-all.php 
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to form and migrate inward to the "standard" location of the X-point (3 — lOMyr). Thus, 
relatively few FU-Ori outbursts will take place after the planets arrive at r = ~ 0.05AU. 

On the other hand, migration will be aided by the fact that the planets might not have 
time to open up gaps in the disk. As a result, we expect that the migration rates will be faster 
than the standard Type II migration rates (for giant planets with gaps). The migration rate 
will be closer to the Type I migration rates, which are fast, and are even faster for larger 
mass planets (that usually open up a gap). 

The second complication is that, because FU-Ori outbursts are short, the star has 
difficulty adjusting its rotation period to co-rotate with the disk: During the quiescent 
phase, the star rotates every ~ 4 days, in sync with the inner disk edge at r ^ 0.05AU. 
During an outburst, the large accretion rate suddenly pushes the disk edge to r ^ 0.02AU. 
For the star to co-rotate with the new disk edge, it needs four times its angular momentum. 
However, the maximum rate for the star to gain angular momentum is J ~ M^rx'^^x- 
For Md = lO-^Moyr-i and rx = 0.02AU, it takes 2000 - 3000 years to spin up the star, 
somewhat longer than the timescale for each outburst. As a result, the inner disk edge 
generally rotates faster than the star, and the magnetic funnel forms a leading spiral. Angular 
momentum flows into the star and makes the truncation radius smaller than its equilibrium 
value (equation [I6]). These departures from equilibrium thus allow migrating planets to 
move even farther inward. 

Next we note that the density of gaseous material in the disk midplane is enhanced 
over the expected value p ~ S/2/7. Because disk truncation by th e funnel flow increases 



the density by a factor r/H ~ 40 at the X-point (jShu et al. 19941 ). the expected density 
px ~ rE/2iJ^. This density enhancement has important ramifications: 

[1] Giant planets that successfully migrate to the X-point at r ^ 0.05AU can be pushed 
further inward (perhaps to r ^ 0.02AU) during FU-Ori outburst events. With the density 
enhancement at the disk midplane, the timescale for inward movement can be short (< 40yr), 
so that migration can take place during a single outburst. 

[2] Planetary cores can migrate inward to the X-point on a short timescale, often fast 
enough that they cannot become gas giants before they reach the inner disk edge. If planetary 
cores are large enough, they continue to accrete gaseous material while they orbit at the X- 
point. The density enhancement considered herein provides a corresponding enhancement 
in the supply of available material. One issue associated with forming hot Jupiters in situ 
is the lack of gaseous material in the inner disk. This gas shortage is largely alleviated for 
planets at the X-point, both from this density enhancement and from the continuing fiow 
of material into the region (Ward 1997b). On the negative side, high temperatures at the 
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X-point make the gas essentially dust-free. Gas accretion onto planets is limited by cooling 
processes, which will be compromised if the opacity is too low. 

4. CONCLUSION 

This letter presents a new migration mechanism for embedded planets. These bodies 
experience a headwind due to the subkeplerian rotation curve of the gas in magnetically 
controlled disks, and drag forces drive inward migration. This mechanism dominates over 
Type I migration (Figure [1]) for sufficiently small planets (mp < IM,^) and/or close orbits 
(r < lAU). The total migration time moderately decreases due to subkeplerian torques, but 
the mass accreted by planetary cores during the migration epoch changes more substantially 
(Figure E]). 

The gap opening process is modified in subkeplerian disks because the planet moves 
faster than the gas (supersonically) both inside and outside its orbit. Gap opening thus 
requires nonlinear interactions (approximately given by equation [10]). The estimated mass 
threshold for gap opening is mp > lOOM^, and the planet resides near the inner gap edge. 

For typical parameters, magnetic truncation of the inner disk occurs at rx ~ 0.05AU, 
and thus naturally explains observed extrasolar planets in ~ 4-day orbits. Migration of 
planets to smaller orbits, as sometimes observed, requires an additional mechanism. During 
FU-Ori outbursts, the X-point moves inward to the required location rx ~ 0.02AU due to 
the large accretion rate, which also allows fast inward migration. X-wind theory (S07,C08) 
thus provides a viable explanation for all of the currently observed short-period planets 
(long-period planets require additional considerations). 

Finally, planetary cores that halt migration at the X-point {rx ~ 0.05AU) can be sub- 
jected to a significantly enhanced gas density and a continuing supply of material (including 
inwardly migrating solids). These features help proto-hot-Jupiters finish forming in their 
observed ~ 4-day orbits. 

This letter shows that magnetic effects drive subkeplerian migration and produce final 
orbits with r ~ 0.02 — 0.05AU (from disk truncation). These processes must compete 
with other mechanisms acting on embedded planets, including Type I migration (Ward 
1997ab), stochastic migration from turbulence (Nelson & Papaloizou 2004, Laughlin et al. 
2004b, Nelson 2005, Johnson et al. 2006, Adams & Bloch 2009), runaway migration (Masset 
& Papaloizou 2003), torques exerted by toroidal fields (Terquem 2003), opacity transitions 
(Menou & Goodman 2004), and radiative (non-barotropic) effects (Paardekooper & Mellema 
2006, Paardekooper & Papaloizou 2008). Future work should study how these processes 
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interact with each other to ultimately form giant planets. 
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